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INTRODUCTION 

The  structural  optimisation  problem  may  be  posed  hy  the  following: 


min  f(x.)  1 

subject  to  the  conditions 


1,2. ..I 


8  (x^J  ^  0  j  "  1,2. ..v 


1  xi  i  xi  (- 

For  the  minimum  weight  design  of  structures  modeled  by  bar  and  membrane 

o 

plate  elements  the  equations  may  be  given  by: 


A.  x. 
x  i 


g.(x  )  -  [  £  E  ,/x.  -  U,  /U,  <_  0. 
J  j=,]_  xj  1  3  3 


Resizing  i9  generally  accomplished  in  Mathematical  Programming  pro¬ 


cedures  by  letting 


xf 1  -  x*  +  nr  S*  (6) 

where  ar  the  step  size  is  usually  selected  arbitrarily  and  s£  is  the 
direction  of  movement. 

Optimality  criteria  (OC)  methods  resize  the  structure  based  on  a  solution 
of  the  problem 

+  E  X  g  -  0  (7) 

by  methods  of  successive  substitutions  where 

x^4"1  -  x*  (n  Xj)  jeJA.  (8) 


Here  n.  is  "resizing  parameter"  and  is  a  set  of  "active"  constraints. 

For  this  report  the  meaning  of  the  n  parameter  is  that  of  the  similar 
symbol  of  Ref.  5.  Procedures  which  compute  the  values  of  the  X  set  con¬ 
sidering  the  equations  for  X  as  coupled  will  be  referred  to  as  "generalized" 
OC  methods'1’  and  those  where  the  X  are  assumed  uncoupled  as  "simple"  OC 


methods  . 


II.  OBJECTIVES 


There  are  two  distinct  projects  associated  with  this  research.  These 


1.  A  preliminary  investigation  of  the  effectiveness  of  a  strategy 
to  reduce  divergence  and  oscillation  in  large  scale  optimal 
structural  synthesis  procedures. 


are: 


2.  The  development  of  a  primal  Mathematical  Programming  (MP)  pro¬ 
cedure  suitable  for  large  scale  structural  synthesis . 

III.  OSCILLATION  AND  DIVERGENCE  CONTROL 

111.1  Background.  Almost  all  optimization  methods  suitable  for  large  scale 

structural  synthesis  require  the  selection  of  a  constraint  set  for  inclusion 

into  the  resizing  problem.  This  is  usually  done  by  including  all  those  con- 
2  3  5 

s train ts  where  *  ’ 

8j  i  .  (9) 

Here  e.  will  be  referred  to  as  the  constraint  "band  width".  The  inclusion 
J 

of  too  many  constraints  results  in  an  unnecessary  increase  in  the  compu¬ 
tational  effort  required  for  the  solution  of  the  resizing  problem.  For 
example,  this  effort  can  be  substantial  and  may  greatly  exceed  the  reanalysis 
effort  in  general  OC  procedures^.  Excess  constraints  may  also  "over cons  train" 
the  problem  producing  heavier  designs. 

On  the  otherhand  selection  of  band  widths  that  are  too  small  may  lead 

to  the  major  violation  of  a  constraint  that  was  not  included  in  the  resizing 

2 

problem  resulting  in  an  increase  in  weight  after  scaling  or  other  boundary 
restoration. 

A  major  difficulty  in  the  use  of  optimization  methods  is  that  no 
rigorous,  or  even  reasonably  reliable,  efficient  procedure  has  been  formu¬ 
lated  for  band  width  selection.  The  desire  to  reduce  the  computational  effort 
and  avoid  the  overconstrai.it  problem  usually  results  in  selection  of  narrow 

band  widths  and  thus  in  occasional,  or  even  frequent,  problems  with  os¬ 
cillation  or  divergence^ ’  . 

Oscillation  or  divergence  resulting  from  such  problems  will  be  referred 
to  here  as  "primary  oscillation".  In  addition  to  this  mode,  oscillation 
can  result  from  too  large  a  value  of  the  resizing  parameter  or  step  size. 

Such  oscillation  will  be  referred  to  as  "secondary  oscillation". 

111. 2  Procedure  for  dealing  with  primary  oscillation.  The  basic  concept 


here  is  to  introduce  a  potentially  active  constraint  set  where  jeJ^  if 

-lj  'e2j  aW 

*2j  '  "lj  (U) 

These  constraints  are  not  included  directly  in  the  resizing  problem. 

However,  if  after  resizing  it  appears,  based  on  the  gradient  information 

associated  with  these  constraints,  that  any  of  them  may  be  violated  the 

resizing  step  is  shortened  in  an  attempt  to  avoid  this  violation. 


Where 


(12) 


4x. 


i+l  r 
xi  ~  xi 


rfl 


the  estimated  value  of  g  ' *  Is  given  by 

t+1  r  .  _  Jr  .  r 
84  *  84  +  ^  8j 


If  any 


5j 


t+1  „ 

e  >0 

8j 


9j  »i  AXi 


compute  a  step  shortening  quantity  K  such  that  a  move  would  produce 


value  of  g 


r+1 

i 


0  for  these  constraints.  Thus 
r 


Kj  "  8j,i  Xi 


Call  the  smallest  of  these  K  and  redefine  as 

xf1  -  K*  (x*  +  xj) 


(13) 

(14) 


(15) 

(16) 


III. 3  Procedure  for  dealing  with  secondary  oscillation.  If  a  weight  in- 

crease  results  after  resizing  and  the  sets  and  Jp  contain  no  con- 

r  r 

straints  that  were  not  in  either  JA  or  Jp  then  it  is  assumed  that  the 
resizing  parameter  is  too  large.  Design  r  +  1  then  discarded  and  a  new 
resizing  move  made  where 
n 

or 


r*'1-  r/2. 


(17) 


am  -  r/2  . 

III. 4  Termination  Procedure.  The  resizing  process  is  terminated  if 
r+1 


IfCx1*1)  -  f(xr)[/f(xr  <  Cx 


(18) 


where  is  a  convergence  criteria  or  if 

n  <  C2  (19) 

where  is  the  minimum  resizing  parameter  size. 

III. 5  Added  Storage  Requirements  and  Computational  Effort.  Once  a  step 
reduction  value  is  computed  all  the  gradient  information  associated 
with  that  constraint  may  be  discarded.  Thus  it  is  only  necessary  to  store 
one  set  of  constraint  gradient  components  rather  than  the  gradients  of 
all  constraints  in  the  potentially  active  set.  Furthermore,  in  many 
procedures  once  the  quantities  associated  the  solution  of  the  resizing 
problem  are  computed  the  gradients  of  the  active  set  are  no  longer  needed. 
Some  part  of  the  active  constraint  gradient  storage  may  therefore  be  used 
for  the  particular  potentially  active  constraint  being  considered.  Thus, 
no  additional  storage  is  required  for  the  primary  oscillation  control  pro¬ 
cedure  in  such  instances. 


In  Che  case  of  OC  methods  using  the  secondary  oscillation  control  pro¬ 
cedure  if  one  is  to  formulate,  a  new  resizing  problem  at  a  point  x^  using 
a  smaller  n  computational  effort  may  be  reduced  if  information  associated  with 
the  1  set  of  the  former  resizing  problem  is  saved.  This  information  is 
derived  from  the  gradients  of  the  active  constraint  set.  Since  active 
constraint  set  gradient  storage  is  no  longer  needed  this  resizing  problem 
information  may  be  stored  in  its  place.  Furthermore,  since  the  active 
constraint  set  gradients  require  much  more  storage  than  the  resizing  infor¬ 
mation  this  extra  storage  may  be  used  for  some  or  all  of  the  potentially 
active  constraint  gradients.  Thus,  with  appropriate  information  storage 
these  procedures  should  normally  not  normally  require  a  substantial  increase 
in  storage  capacity. 

The  calculation  of  the  additional  potentially  active  constraint  gradients 
requires  additional  computational  effort  at  each  redesign  cycle.  This 
additional  effort  is  however  usually  small  compared  to  the  total  compu¬ 
tational  effort  and  much  smaller  than  the  effort  wasted  by  bad  moves  re¬ 
sulting  from  failure  to  consider,  at  all,  these  potentially  active  con¬ 
straints. 

IV.  MATHEMATICAL  PROGRAMMING  (MP)  ALGORITHM 
IV. 1  Background.  Mathematical  programming  procedures  are  considered  to  be 
the  most  general  and  easily  applied  of  the  various  optimization  procedures. 

g 

They  avoid  some  of  the  difficulties  associated  with  OC  procedures  .  MP 

procedures  are  however  generally  considered  to  be  poorly  suited  to  large 

structural  synthesis  problems  without  use  of  approximation  techniques, 

at  least  in  their  primal  form,  due  to  the  relatively  large  number  of 

4 

reanalysis  typically  required  for  most  procedures  . 

It  has  been  the  feeling  of  this  writer  that  the  view  that  MP  procedures 
are  too  inefficient  for  large  scale  structural  synthesis  has  some  justifi¬ 
cation  in  the  case  of  most  existing  methods.  It  was  also  felt,  however, 

that  should  be  relatively  easy  to  construct  a  procedure  based  on  the  method 

q 

of  feasible  directions  which  will,  on  the  basis  of  total  computational 
effort,  be  competative  with  OC  procedures.  Such  a  procedure  would  provide 

g 

the  flexibility  to  treat  problems  poorly  suited  to  OC  methods  .  An  attempt 
at  such  a  procedure  is  presented  here. 


IV. 2  The  Feasible  Direction  (FD)  MP  Procedure.  The  feasible  direction 
problem  is  usually  formulated:^"  Find  and  a  so  as  to 

Maximize  a  (20) 


where 


SAVf  Cxi)  +  a  <_  0 


SX7g  0c±)  +  W  a  <  0 


-1  <  S.  <  1  (23) 

where  S  consiting  of  components  is  the  "feasible  direction"  of  move¬ 
ment,  a  is  a  dummy  variable  and  a  weighting  parameter. 

Equations  (20)  and  (21)  state  that  on  the  basis  of  the  linearized 
functions  the  solution  to  this  problem  will  produce  a  maximum  possible  im¬ 
provement  in  f(x^).  Equations  (20)  and  (22)  state  for  W.  =  0  there  will 
be  no  movement  toward  constraint  violation  and  where  >  0  there  will  be 
movement  away  from  violation.  Equation  (23)  is  used  to  eliminate  unbounded 
solutions.  This  Is  a  linear  programming  problem  and  may  efficiently  be 
solved  by  one  of  many  well  developed  procedures. 

Here  one  has  the  problem  of  determining  an  appropriate  active  constraint 
set  to  include  in  Eqn.  (22) .  Too  large  a  band  width  will  overconstrain  the 
problem  by  forcing  movement  essentially  parallel  to  or  away  from,  a  con¬ 
straint  that  may  not  be  critical.  Too  small  a  band  width  may  produce 
serious  violation  of  a  constraint  that  was  not  considered  in  the  direction 
finding  problem  of  Eqns.  (20-23). 

The  weighting  parameter  in  most  feasible  directions  methods  is  given 
as  positive  in  order  to  avoid  constraint  violation.  Such  violation  will 
occur  on  the  convex  constraints  which  are  usually  encountered  in  structural 
design.  The  central  idea  of  the  improved  method  is  to  set  such  that  it 
will  produce  movement  toward  the  constraint  if  one  is  near  but  not  on  a 
critical  constraint.  This  eliminates  the  overconstraining  effect  of  too  large 
a  band  width.  Thus  replace  the  a  term  in  (22)  with  a  term  producing 
movement  toward  the  constraint  and  thereby  rewrite  Eq  (22)  as 

«r  STVgJ  (x±)  <  gj  .  (24) 

It  may  be  seen  that  left  hand  side  will  produce  an  estimated  change  at 
most  equal  to  the  constraint  value.  Then  after  resizing  the  value  of  the 
constraint  will  move  toward  zero  or  even  near  zero  if  the  constraint 


is  critical.  In  otherwords,  if  movement  toward  violation  of  the  constraint 
improves  the  design,  Eq.  (24)  will  limit  that  movement  only  to  the  extent 
required  to  avoid  violating  the  constraint.  Thus,  the  redesign  move 
will  be  influenced  only  by  constraints  expected  to  be  critical. 

It  should  Be  noted  this  procedure  may  also  produce  oscillation  or 
divergence  resulting  from  deficiencies  in  band  width  or  n  selection  the 
same  procedures  described  earlier  for  control  of  oscillation  and  divergence 
may  of  course  be  used  here. 

The  value  for  ar  for  this  study  is  given  by 

ar  -  nr|sr"1iVfr  (25) 


where 


I  and  n  is  arbitrarily  selected.  Such  a  value  of  h  will 


tend  to  produce  a  change  in  the  objective  function  of  approximately  100b 
percent  if  an  S  move,  similar  to  S  ,  is  made  in  the  Vf  direction.  Thus 
n  »  .5  would  tend  to  produce  about  a  50%  change  in  weight  if  one  moved 
in  the  Vf  direction.  Since  the  direction  finding  problem  would  produce 
deflected  away  from  Vf  by  the  constraints,  the  change  after  an  actual  move 
would  be  substantially  less  than  that  produced  by  a  move  in  the  Vf  direction. 
This  would  be  particularly  true  during  the  later  stages  of  the  search  were 
the  design  is  highly  constrained. 

If  a  redesign  move  fails  to  produce  a  weight  reduction  n  is  reduced 
and  design  x  is  discarded  as  in  Section  III. 3.  A  new  design  is  then 
generated  using 

n**1  -  nr/2  .  (26) 

When  this  occurs  the  band  widths  are  also  narrowed  by  setting 

ef*"1  -  eJ/2  (27) 

o  3  3 

The  e  are  arbitrarily  selected.  Similarly  these  cuts  are  made  if  S  ■  0 

3  9 

is  the  solution  . 

The  design  is  restored  to  the  feasible-infeasible  boundary  by  scaling 

as  in  Ref.  2  or  by  some  other  procedure  where  the  problem  is  not  of  the  form 

9 

of  eqns  (4-5)  . 

Termination  is  accomplished  by  the  procedure  of  Section  III. 4. 

V.  EXAMPLES 

2 

V, 1  Problems.  The  two  ten  bar  truss  problems  posed  by  Venkayya  et  al 
are  repeated  in  this  preliminary  study  using  the  parameters  of  Ref.  5. 

The  stress  "constrained”  problem  involves  the  minimum  weight  design  of  the 


an  indeterminate  structure  under  a  single  loading  condition  where  the  stress 
in  all  members  must  be  held  at  or  below  one  of  two  specified  values.  The 
optimal  design  in  this  structure  weighs  1,497.60  lbs.  In  the  displacement 
constrained  problem  two  additional  constraints  are  placed  on  the  deflection 
of  two  joints  and  only  one  stress  limit  value  is  used.  Side  constraints 
[eqn  (3)]  are  used  for  all  variables. 

The  stresses  and  deflections  are  determined  by  a  finite  element 

analysis  the  results  of  which  may  be  used  to  develop  the  constraint  gradient 

2 

information  by  means  of  the  virtual  unit  load  method  . 

These  are  now  classic  benchwork  problems  and  represent  a  moderately 
difficult  challenge  for  any  proposed  structural  optimization  scheme.  The 
stress  constrained  problem  has  8  (of  10)  active  stress  constraints  in 
addition  to  several  active  side  constraints.  The  displacement  constrained 
problem  has  two  local  minima  one  of  which  is  constrained  by  the  two  dis¬ 
placement  constraints  (this  design  weighs  5,077.6  lbs)  and  one  by  a  stress 
and  one  displacement  constraint  (this  design  weighs  5,061  lbs).  Several 
side  constraints  are  critical  in  both  local  optima. 

V. 2.  Procedures.  The  generalized  OC  procedures  described  by  Khot  et  al~* 
as  the  exponential  [eq  (7)  of  Ref.  5],  and  linear  [eq  (9)  of  Ref.  5]  re¬ 
cursion  forms  was  modified  to  use  the  oscillation  control  techniques 
described  here  in.  In  addition  this  procedure  was  also  adapted  to  use  the 
multiple  iteration  Newton  Raphson  procedure  described  by  Khot  et  al  in 
Ref.  7.  The  procedure  from  Ref.  7  is  remarkably  similar  in  that  of  Ref.  4 
except  for  the  method  of  solution  of  the  \  problem.  The  oscillation  control 
methods  were  added  to  the  experimental  program  used  to  obtain  the  results 
given  in  Ref.  5  and  7  to  generate  the  results  contained  herein. 

The  MP  program  used  here  was  also  developed  from  this  same  experiment¬ 
al  program  by  replacing  the  resizing  portions  of  the  program.  It  did  not 
employ  the  primary  oscillation  control  technique. 

71.  RESULTS  OF  USE  OF  OSCILLATION  AND  DIVERGENCE  CONTROL  PROCEDURE 

VI.  1  Primary  Oscillation  Control.  Table  1  illustrates  the  application  of 
this  procedure  to  two  cases  where  serious  oscillation  was  experienced 
using  OC  procedures^ *7.  Both  involve  the  stress  constrained  problem.  In 
the  first  case  the  exponential  recursion  form  with  multiple  iterations 
solution  for  the  X  set  was  employed  and  in  the  other  the  exponential  form 


with  single  solution  for  the  set  is  used.  It  may  be  seen  that  the  use 
of  this  procedure  did  in  fact  control  primary  oscillation. 

Two  other  test  cases  were  used  in  which  OC  the  procedures  described 
in  Ref.  5  and  7  were  induced  to  oscillate  badly  by  use  of  an  excessively 
large  resizing  parameter  (n  a  1) .  The  use  of  the  primary  oscillation  con¬ 
trol  again  suppressed  this  oscillation  mode. 

VI. 2  Secondary  Oscillation  Control.  Table  2  illustrates  the  application 
of  the  secondary  procedure.  With  only  primary  control  one  sees  here 
secondary  oscillation  at  design  number  4-6  and  divergence  at  design  number 
20.  It  may  be  seen  that  the  reduction  in  the  resizing  parameter  produces 
convergence  to  the  optimal  design. 

This  procedure  was  found  satisfactory  only  with  the  linear  recursion 
forms.  It  was  found  unsatisfactory  for  the  exponential  form  since  the 
assumption  that  the  X  problem  formulated  on  the  basis  of  the  linear  re¬ 
cursion  relation  may  be  used  with  the  exponential  form  was  not  valid  at 
small  n  values.  Because  of  this  all  subsequent  numerical  experiments 
are  based  on  the  linear  recursion  forms. 

All  runs  used  e^  *  2  e^  to  define  the  potentially  active  constraints. 

The  quantities  e7 ,  for  the  active  constraints  are  defined  by  the  procedure 

_6 

of  Ref.  5.  Termination  constraints  used  for  all  runs  were  =*  10  ,  *  0. 

VI. 3  Use  of  a  Large  Initial  Resizing  Parameter.  Since  the  early  experiments 
indicated  that  these  controls  would  inhibit  oscillation  that  would  usually 
occur  with  excessive  value  of  n  an  experiment  was  performed  to  investigate 
the  possibility  of  using  a  large  initial  value  of  this  parameter  so  as  to 
speed  convergence.  The  result  is  shown  in  Tables  3  and  4.  These  results 
fail  to  support  the  hypothosis  that  a  large  initial  n  is  advantageous. 

VI. 4  Convergence  of  Multiple  and  Single  Problem  Solutions  Procedures.  It 
may  also  be  seen  from  Tables  3  and  4  that  the  multiple  X  iteration  procedure 
possesses  much  better  convergence  properties  than  the  single  X  solution  pro¬ 
cedure.  In  fact  the  single  X  solution  procedure  fails  on  the  displacement 
constrained  problem.  It  also  fails  on  this  problem  when  no  oscillation 
controls  are  used. 

An  oscillation  control  procedure  is  important  in  allowing  exploita¬ 
tion  of  the  multiple  X  iteration  approach  because  of  the  oscillation 
problems  with  this  method  encountered  in  earlier  studies.  It  is  interesting 


TABLE  2.  Design  Sequence,  Stress  Constrained  10  Bar  Truss  Problem  Single 
Solution,  Linear  Form  n  ■  0,5  at  Start,  Weight  in  lbs. 


Design  No. 

Primary  Control  Only 

Both  Controls 

1 

3,435 

3,435 

2 

2,444 

2,444 

3 

1,836 

1,836 

4 

2,294 

2,294* 

5 

1,863 

1,726 

6 

2,061 

1,653 

7 

1,727 

1,586 

8 

1,718 

1,556 

9 

1,569 

1,539 

10 

1,541 

1,529 

11 

1,524 

1,523 

12 

1,518 

1,518 

13 

1,513 

1,513 

14 

1,509 

1,509 

15 

1,506 

1,506 

16 

1,503 

1,503 

17 

1,501 

1,501 

18 

1,499 

1,499 

19 

1,497.6 

1,497.74 

20 

1,656 

1,656* 

21 

1,685 

1,572* 

22 

1,681 

1,534* 

23 

1,677 

1,516 

24 

1,677 

1,497.67 

25 

1,678 

1,506* 

26 

1,497.63 

27 

1,502* 

28 

1,497.61* 

♦Halved  resizing  parameter  at  this  redesign  cycle 


TABLE  3. 

Comparison  of  Design  Sequences  Using  Different  Initial  Resizing 
Parameters,  Stress  Constraints  10  Bar  Truss,  Weight  in  lbs. 

Design  No. 

Multiple  X  Iterations 
rt*l  at  Start  n*0. 

.5  at  Start 

Single  X 
n-1  at  Start 

Solution 
h*0.5  at  Start 

1 

3,435 

3,435 

3,435 

3.435 

2 

2,128 

2,059 

2,740 

2,444 

3 

1,915 

1,959 

1,825 

1,836 

4 

1,772 

1,801 

1,719 

2.294* 

5 

1,638 

1,664 

1,628 

1,726 

6 

1,593 

1,654 

1,580 

1,653 

7 

1,532 

1,540 

1,548 

1,586 

8 

1,527 

1,532 

1,531 

1,555 

9 

1,521 

1,526 

1,519 

1,539 

10 

1,512 

1,521 

1,509 

1,529 

11 

1,508 

1,516 

1,499 

1,523 

12 

1,505 

1,512 

1,517* 

1,518 

13 

1,502 

1,508 

1,638* 

1,513 

14 

1,500 

1,505 

1,497.9 

1,509 

15 

1,558* 

1,502 

1,572* 

1,506 

16 

1,497.62 

1,499 

1,502 

1,503 

17 

1,497.60 

1,498 

1,497.7 

1,501 

18 

1,497.60+ 

1,532* 

,  1,516* 

1,499 

19 

1,497.63 

1,506** 

1,498 

20 

1,497.60 

1,656* 

21 

1*497.60+ 

1,572* 

22 

1,534* 

23 

1,516* 

24 

1,497.67 

25 

1,506* 

(28)1,497.61+ 

*Halved  resizing  parameter. 
tTerminated  by  convergence  specification. 

♦★Terminated  by  minimum  resizing  parameter  specification 


TABLE  4.  Comparison  of  Design  Sequences  Using  Different  Initial  Resizing 
Parameters,  Displacement  Constrained  10  Bar  Truss,  Weight  in  lbs 


Design  No. 

Multiple  A  Iterations 
n*l  at  Start  '  n=Q.5  at  Start 

Single  A  Solution 
nal  at  Start  r\“0.5  at  Start 

1 

8,266 

8,266 

8,266 

8,266 

2 

6,893 

6,646 

7,146 

6,667 

3 

5,816 

5,824 

6,932 

6,196 

4 

5,593 

5,703 

6,641 

5,947 

5 

5,466 

5,597 

5,376 

5,759 

6 

5,327 

5,471 

5,221 

5,628 

7 

5,358 

5,353 

5,216 

5,516 

8 

5,505* 

5,195 

5,162 

5,401 

9 

5,200 

5,079.6 

5,179 

5,280 

10 

5,094 

5,078.9 

5,129 

5,187 

11 

5,076 

5,084.4* 

5,118 

5,143 

12 

5,066 

5,077.0 

5,115 

5,115 

13 

5,061.0 

5,076.7 

5,113 

5,108 

14 

5,060.87 

5,076.7 

5,112 

5,090.4 

15 

5,060.85 

5,117 

5,104.4* 

16 

5,060.85 

5,111.5 

5,094.2 

17 

5,111.35 

5,089.4 

18 

5,111.31 

5,088.0 

19 

5,111.253 

5,086.9 

20 

5,111.247 

5,086.0 

25 

5,098.7* 

5,083.6* 

30 

5,094.6 

5,082.9+ 

35 

5,093.4 

50 

5,093.3 

*Halved  resizing  parameter. 
tTerminated  by  convergence  specification 


TABLE  5 


Comparison  of  Inverse  and  Ordinary  Variables  Using  Multiple 
Iteration  Linear  Form  Procedure,  n  -  0.5  at  Start,  Weight  in  lbs 
10  Bar  Truss 


Desien  No. 

Inverse 

Ordinary 

Inverse 

Ordinary 

1 

3.435 

3,435 

8,266 

8,266 

2 

2,098 

2,059 

6,912 

6,645 

3 

1,903 

1,959 

6,438 

5,824 

4 

1,662 

1,801 

6,194 

3,703 

5 

1,624 

1,664 

6,019 

5,597 

6 

1,537 

1,654 

5,877 

5,471 

7 

1,518 

1,540 

5,750 

5,353 

8 

1,513 

1,532 

5,628 

5,195 

9 

1,509 

1,526 

5,511 

5,080 

10 

1,503 

1,521 

5,402 

5,078.9 

11 

1,501 

1,516 

5,307 

5,084* 

12 

1,593* 

1,512 

5,182 

5,077.0 

13 

1,497.67 

1,508 

5,101 

5,076.7 

14 

1,497.60 

1,505 

5,094 

5,076.7 

15 

1,565* 

1,502 

5,081 

16 

1,497.60 

1,498.8 

5,094 

17 

1,497.9 

5,081 

18 

1,532* 

5,0770 

19 

1,497.63 

5,076.7 

20 

1,497.60 

5,107* 

21 

1,497.60 

5,076.7 

♦Halved  resizing  parameter. 
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VII.  PERFORMANCE  OF  THE  FEASIBLE  DIRECT I ONS  ALGORITHM 


The  results  of  this  algorithm  on  the  example  problems  using  various 
values  of  n°  are  shown  in  Table  6.  In  all  cases  e°  =  0.5.  A  small  n° 
produces  small  changes  in  the  initial  designs.  A  large  n°  on  the  other- 
hand  produces  early  oscillation  associated  with  the  need  for  step  size 
reduction.  This  of  course  is  similar  to  the  situation  in  0C  procedures. 

Further  development  can  undoubtably  substantially  improve  performance. 
For  example  the  algorithm  described  in  section  V  was  modified  such  that; 

1)  the  equations  (21)  and  (28)  are  invoked  (step  size  is  reduced)  after 
the  change  in  weight  on  scaling  is  greater  than  the  net  decrease  in  weight 
after  redesign  and  scaling;  2)  the  S  in  eq.  (25)  replaced  by  S*  where 
S*  is  obtained  after  iterative  solution  of  the  feasible  direction  problem 
[eqns.  (19-25)]  until  ] S* | = | Sr | .  This  modification  produced  a  dramatic 
improvement  in  performance  (see  Table  6)  on  the  sole  problem  on  which  it 
was  tes  ted . 

VIII.  COMPARISON  OF  THE  PC  AND  FD  PROCEDURES 

On  the  basis  of  this  limited  study  if  one  uses  similar  values  of  the 
resizing  parameter  it  appears  that  the  MP  procedure  requires  fewer 
reanalyses  for  convergence  than  the  single  X  solution  OC  procedure  even 
after  the  oscillation  control  improvements  are  made  to  the  latter.  The 
multiple  A  iteration  approach  seems  to  possess  superior  convergence  prop¬ 
erties  when  these  controls  are  used.  It  appears,  however,  that  it  may  be 
relatively  simple  to  greatly  improve  the  performance  of  the  MP  procedure 
to  the  point  where  it  is  comparable  to  the  multiple  X  iteration  OC  pro¬ 
cedure  . 

On  the  basis  of  computational  effort  required  for  convergence  the 
picture  is  somewhat  different.  Analaysis  and  resizing  times  are  shown 
in  Table  7.  The  analysis  time  is  the  CPU  time  required  to  do  the  finite 
element  analysis.  The  resizing  time  includes  the  time  required  to  compute 
the  necessary  derivatives  and  set  up  and  solve  either  the  X  or  the  feasible 
direction  problem. 

It  may  be  seen  that  the  resizing  time  required  for  the  MP  procedure 
is  similar  to  the  single  X  solution  OC  method  which  is  much  less  than  the 
multiple  X  iteration  approach.  Furthermore,  it  is  the  resizing  effort 
that  dominates.  On  larger  problems  one  would  expect  a  similar  situation. 


TABLE  6.  Design  Sequences  for  Feasible  Directions  Algorithm  Using  Different 
Resizing  Parameters,  Weight  in  lbs.  10  Bar  Truss. 


Design  No. 

Stress 

0.1 

Constrained 

0.2  0.5 

Displacement  Constrained 

0.1  0.2  0.5  0.5 

1 

3,435 

3,435 

3,435 

8,266 

8,266 

8,266 

8,266 

2 

3,166 

2,885 

2,224 

7,770 

7,394 

7,006 

7,006 

3 

2,884 

2,455 

2,722 

7,040 

6,837 

7,642 

6,610 

4 

2,641 

2,096 

1,818 

6,801 

6,629 

7,608 

6,385 

5 

2,448 

1,817 

2,300 

6,615 

6,402 

6,331 

6,171 

6 

2,279 

2,006 

1,787 

6,466 

6,097 

6,217 

5,976 

7 

2,125 

1,703 

1,677 

6,309 

5,826 

5,796 

5,854 

8 

1,854 

1,693 

1,604 

6,150 

5,668 

5,492 

5,564 

9 

1,736 

1,623 

1,557 

6,002 

6,358 

5,225 

5,398 

10 

1,792 

1,571 

1,528 

5,827 

5,512 

5,164 

5,176 

11 

1,671 

1,535 

1,507.6 

5,637 

5,172 

5,281 

5,115 

12 

1,644 

1,516 

1,501.4 

5,396 

5,114 

6,225 

5,739 

13 

1,624 

1,742 

1,737 

5,178 

5,090 

5,175 

5,066.9 

14 

1,606 

1,506 

1,498.3 

5,112 

5,088 

5,125 

5,126 

15 

1,585 

1,737 

1,498.4 

5,169 

5,090 

5,108 

5,091 

16 

1,567 

1,501 

1,497.9 

5,082 

5,081 

5,099 

5, 06  7. '6 

17 

1,549 

1,725 

1,498.5 

5,079 

5,084 

5,095 

18 

1,524 

1,498 

1,497.7 

5,077 

5,079.6 

5,088 

19 

1,515 

1,725 

1,498.4 

5,078 

5,079.6 

5,085 

20 

1,505 

1,549 

1,497.6 

5,077.0 

5,079.4 

5,081 

25 

1,506 

1,497.6 

5,077.6 

5,078.6 

5,077. 

1 

30 

1,497.6 

5,076.7 

5,077.9 

5,076. 

7 

35  5,077.6 

50  5,076.7 


tModified  Algorithm 


Ten  Bar  Truss 

Stress  Constrained  *022 


1.996 


.022  .1330 


.073  -274 


Ten  Bar  Truss 

Displacement 

Constrained 


.021 


.379  .022  .071 


.022  .112 


Two  Hundred  Bar 


7AQ 


Furthermore,  on  large  problems  a  more  efficient  but  more  complex  linear 
programming  procedure  which  ignores  zero  matrix  entries  rather  than  the 
simple  procedure  used  here  should  substantially  reduce  the  resizing  effort 
making  the  MP  procedure  more  attractive. 

IX.  CONCLUSION 

Much  more  work  needs  to  be  done  to  verify  the  preliminary  results 
developed  here  after  further  refinement  of  the  concepts  presented.  The 
results  of  this  work,  however,  supports  the  initial  assumptions  that  os¬ 
cillation  problems  associated  with  many  optimization  methods  may  easily 
be  greatly  reduced  and  that  a  simple  primal  MP  procedure  without  approxi¬ 
mations  can  be  competitive  with  OC  procedures  for  finite  element  based 
structural  synthesis. 

On  the  basis  of  this  preliminary  study  the  MP  procedure  seems  more 
attractive  than  the  OC  procedures  on  large  problems  with  many  active  con¬ 
straints  due  to  the  large  effort  required  to  set  up  the  X  problem^.  On 
problems  with  very  few  active  constraints  the  ability  of  the  OC  procedures 
to  produce  very  large  initial  weight  reductions  makes  these  procedures 
attractive.  Additional  work  needs  to  be  done  however  to  reduce  the  resizing 
computation  effort  of  the  multiple  X  iteration  approach  to  allow  exploita¬ 
tion  of  its  superior  convergence  properties . 

These  conclusions  are  of  course  only  tentative  and  are  based  on  very 
limited  evidence.  The^  importance  of  these  initial  successes  however  justifies 
expanded  research  on  these  techniques . 

X.  RECOMMENDATIONS 

The  preliminary  results  of  this  early  research  are  quite  encouraging 
and  justify  further  study  since  these  methods,  if  successful,  represent 
major  advances.  The  following  additional  research  is  therefore  reconmended: 

1.  Application  of  the  oscillation  control  techniques  to  a  simple 
OC  procedure. 

2.  Refinement  of  the  techniques.  For  example  the  development  or 
improvement  in  the  methods  of  band  width  and  step  size  deter¬ 
mination  or  specification. 

3.  Treatment  of  additional  static  problem  examples. 

4.  Treatment  of  dynamic  problem  examples. 

5.  Treatment  of  example  problems  with  local  buckling  constraints. 


Comparison  of  the  process  developed  as  the  result  of  this 
research  with  important  large  scale  synthesis  capabilities  sac 
as  tfie  ACCESS3  and  OPTSTAT  codes. 

u  justified,  the  incorporation  of  successful  new  methods  into 
e  formal  structural  synthesis  program  for  general  distribution. 
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